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I. SUMMARY AND INTRODUCTION 



Daring the past decade aircraft speeds have extended further and 

further into the supersonic regime - . Today , most aircraft design 

companies project civil transports to fly between Mach 2 and Mach 3. 

Such machines will experience structural problems due to the aerodynamic 

heating associated with these high speeds. Therefore there is a need 

for structural research into the effects of this type of environment. 

Early studies on radiant heat techniques for thermal simulation were 

[l 2 3 4.5]* 

made by Abramson in the United States and Horton > > > > j. n great 

Britain concurrently. Their works proposed computer controlled systems 

using quartz tubular heaters as the basic heat producing elements. The 

concept of a time-shared computer for the calculation of desired radiant 

[5] 

heat flow was described by Horton . The idea has been amplified and 
enlarged in this work. 



A practical aerodynamic heating computer-controller must be: 

(1) Capable of producing the correct temperatures over a given 
range of Mach numbers. 

(2) Relatively simple to operate. 

(3) Dependable. 

(4) Safe to operate. 

,(5) As inexpensive as possible. 



To create the thermal environment corresponding to a given Mach 
number, the simulator must produce the appropriate heat flow. This has 
three components. 

(1) The convective heat transfer to the body. 

(2) The radiant heat transfer to or from the body. 

(3) The internal heat generated within the body. 



■X* 

Superscript numerals in brackets refer to the References at the 
end of the report. 
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It is normally accepted that the Newtonian law can be used to compute 
the convective heat flow. Thus, the convective heat input condition 
is defined by 

q = h(t - t ) 
u aw s 

where :• 

2 

q = heat flow per unit area (BTU/ft ) 
h = heat transfer coefficient (BTU/hr/ft^/— F) 
t = adiabatic wall temperature (°F abs) 

t = suface temperature (°F abs) 

The purpose of this work was to design and build a control system 
to meet the above conditions. Even though the equipment is the result 
of a university research project, it is felt that many of the ideas 
developed are directly applicable to industrial plants. 

The results of the investigation show that the use of a time-shared, 
digital-analog computer control is feasible and most desirable from a 
standpoint of cost, accuracy, simplicity and dependability. Although 
the control was limited to a four channel heating system, the results 
show that it is possible to expand the system to more channels, and to 
include loading should this be desired. 
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II. BASIC LOGIC AND BLOCK DIAGRAM 



As stated previously the accepted Newtonian equation: 

q = h(t - t ) 
aw s 



was used to determine the desired convective heat flow, q , into the 

body. External radiant heat input, q^ (such as from a bomb flash or 

solar radiation) was considered to be a quantity which could be positively 

defined and supplied as an analog from a function generator. (in fact, 

the basic digital to analog converter was adapted for this purpose. ) 

Radiant heat loss in flight was considered as a correction term to the 

4 

q supplied, since the computation K n t and the computation 

J- s 



dt 



K, 



2 citT 



are closely allied. 

The quantities t Qw , h , and q^ can, for any specific flight 
program, be defined explicitly as time dependent variables. Therefore 
analogous time dependent voltages can be produced from appropriate func- 
tion generators. However, it is our premise that if the operations 
involving t gw , h , and q^ can be sequentially programmed, then a 
single digital to analog converter can be adapted to generate all. 

(See Fig. l) In the system developed, a tape reader *is used to set up 
a voltage analogous to q^ at time t in the program. This voltage is 
stored in a simple capacitor store using a peak follower circuit as 
described later. Next, at the time t + At the adiabatic wall temper- 
ature, t , is generated. The appropriate surface temperature, t 
aw [7] s 

is measured by a thermocouple and is combined with in a summa- 

tion circuit described in Section III. The output from the summer is 
stored and when the tape advances to read the heat transfer coefficient, 
the stored temperature potential is fed into the converter as a reference 
voltage. The converter then multiplies the temperature potential by the 
heat transfer coefficient, thus computing q . At the same time q^ , 
previously stored,, is added to q giving a voltage analogous to the 
desired heat flow, q^ , uncorrected for radiant loss. The signal is 
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fed into a comparator where it is compared with the measured heat flow 
corrected for radiant loss. The difference, C , between the computed 
and measured values is used to adjust the output of the saturable core 
reactors so that C converges to the control point for correct heat 
flow. Details of the comparison system, the saturable core reactor 
control, and the heat input determination are given in Ref. [ 8 ] and [ 9 ] 
respectively. 

The basic logic shown in the block diagram (Fig. l) was used in the 
design of the control system discussed. In order to satisfy the design 
requirements as previously stated, the system incorporates time sharing 
and digital to analog conversion techniques . Figure 2 is a diagram of 
the time-shared computer control system. In order to familiarize the 
reader with its operation, the time-wise sequence of events required to 
compute q^ , will be followed. 

(The system components are described in Section III, and Appendix 
A contains basic analog computer techniques . ) 

First, assume the wipers are on pole one of each of the nine banks 

and the tape input command is q_^ . The reader generates a voltage 

proportional to q^ from the reference voltage and this is fed through 

pole one of bank one to the storage capacitor C . The switch is now 

stepped to position two and the tape command is simultaneously changed 

to t . The reader generates a voltage analogous to t^ from the 

original reference voltage, t , as a voltage, passes through pole 
* aw 

two of bank one to the storage capacitor . The output of peak 

follower two is analogous to t and, via summer four, has t sub- 

aw * s 

tracted. This quantity, -(t - t ), is now an analog signal appearing 

aw s 

at pole three of bank one. Once again the switch advances and the tape 

command changes to h . The analogous voltage to -(t “ t s ) appears 

as a reference voltage and is multiplied by the command signal in the 
reader head giving q . At the same time, q^ is added to q in ampli- 
fier one giving the analogous voltage of q^ . This voltage passes 
through pole three of bank six, seven, and eight/ and is stored in C-^ a 

and . The output of peak follower six is a continuous voltage 

analogous to q^ . As the steps continue-, channels two, three, and 
four are filled in order. 
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Banks four and five are used to empty the storage capacitors in 
order to allow the storage of decreasing voltages. 

Banks seven and eight isolate . This is necessary in order 
to prevent discontinuous signals of q^ and still allow for decreasing 
voltage levels analogous to desired heat flow. 
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III. COMPONENT DESIGN EVALUATION 



An eight pin Soroban tape reader is utilized in the function gen- 
erator section of the control system. The reader is used to connect 
various combinations of resistors to amplifier one. As shown in 
Appendix A, the output voltage of a summer amplifier is : 




Therefore , by proper choice of resistor values, and by correct choice 
of hole combinations in the tape, the reader head in conjunction with 
amplifier one, will produce any multiple of a given voltage. The prin- 
ciple of this operation is illustrated in Fig. 3, a simplified diagram 
of the reader head amplifier combination. The choice of the resistors 
was made to correspond with the values of numbers on a Soroban tape 
puncher used in programming. Since the tape puncher keyboard is limited 
to numbers from 1 to 159.> the maximum inaccuracy of generation is 

, x 

2(159) 

where X is the operating reference voltage range. If full utilization 

of the eight holes had been made in the tape puncher, 255 combinations 

would result, with a maximum inaccuracy of 
< 

± X. 

2(255) 

Therefore, unless conformity to available equipment is a criteria, any 
accuracy desired may be obtained, since the inaccuracy is equal to 

± x 

2(2 n - 1) 

where n is the number of holes. 

Amplifiers used in the digital-analog computer control are Model 
K2-XA operational amplifiers made by Philbrick Researches Incorporated. 
The K2-XA has a small signal rise time of 1 psec , an input impedance 
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of over 100 Meg ft, output range of ± 100 v , and a drift rate referred 
to the input of approximately 8 mv per day. 

Amplifier 3, Fig. 2, is a chopper stabilized amplifier with a DC 
gain of 1000. This amplifier uses a Philbrick K2-P stabilizing amplifier 
in conjunction with a K2-XA in order to increase stability, giving it a 
long term drift rate of less than 0.1 mv/da y referred to the input. 

The peak follower circuit mentioned previously is simply a K2-XA 
amplifier with a feedback loop of zero resistance'. A capacitor, isolated 
from the rest of the circuit by a diode and the amplifier, is charged to 
the peak value of the input signal. The amplifier output produces a con- 
tinuous voltage equal to the charge on the capacitor. Discharge time 
constants of five minutes or longer are obtainable- 

Multi -contact rotary stepping switches were used in the control 
system. Four milliseconds switching time is required with a self- 
excited stepping capability of 58 steps per second. The electrical 
noise resulting from the use of the stepping switch was not a problem 
since signal voltages were kept large whenever possible to completely 
overshadow the microvolt noise level. Also, the reader delays 

inputs until after all wipers have attained their new positions. 

A Philbrick dual power supply, Model R-100B was used to drive all 
the operational amplifiers. Typical drift for this power supply is less 
than 100 mvdc per day. 

«r 

The storage of voltages in the control system is done entirely by 
the use of capacitors. and C shown in Fig. 2 are 0-5df while 

the rest of the storage capacitors are 0.25hf • The charging time of 
the capacitors is ten milliseconds. After this period, the capacitors 
have attained 99-9$ of their final charge. The charging time of 10msec 
would permit switching up to sixty steps per second. The speed in this 
control system was limited to ten steps per second by the motor drive 
on the reader head. Time base for the entire control system was gener- 
ated by an adjustable pulsing circuit shown in Fig. k. This pulse was 
used to synchronously step the reader head and switches. 
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IV. OVERALL SYSTEM EKGIHEERIKG AMD PERFORI-IANCE 



The four channel control system performance data is found in 
Appendix C. Tests were run to determine the accuracy and repeatability 
of multiplication by the computing system. Since accuracy of computation 
is directly effected by the peak follower output; determination of the 
R-C time constants for the storage network was also made. 

Accuracy of multiplication was found to be in error by 1.6$ or less 
for computing speeds of eight steps per second. This accuracy was de- 
termined by multiplying two known reference voltages together; measuring 
the product and comparing it with the otherwise calculated value, (see 
Appendix C) 

Because of the inherent stability of the computer components and 
care taken in wiring to guard against ground loop problems; the repeata- 
bility of computation over a fifteen minute period resulted in output 
variations of less than 1.0$. 

. R-C time constants for the peak follower inputs were found to be 
five minutes or greater resulting in peak follower output variations of 
less than 0-5$ for a switching time of 0.12 seconds. 

Accuracy of producing q , h ; and t is dependent upon the 
number of discrete levels of the reader; the tolerances of the resistors 
used; and the accuracy of the reference voltage. First, 159 numbers were 
available giving the reader a maximum inaccuracy of 0 . 3 $. The reference 
voltage was a dry cell, and since the current demand was at most a few 
microamperes, it was considered to have constant output. Finally, the 
measured error due to resistor tolerance was ± 3 $ maximum. 

As stated previously, maximum switching speed was limited to 10 
cycles per second by the motor drive on the reader head- Since three 
steps are necessary for each determination of q^ , corrected outputs 
of desired heat flow were made every 1-2 seconds for each of the four 
channels. This speed could have been increased by a factor of five 
simply by changing the motor speed without adversely effecting the ac- 
curacy of the control system. 
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Total cost of the control system described in this work was 
approximately $ 1 , 500 - 

To prevent the hazards of infra red radiation exposure to operating 
personnel, the control system and instrumentation was remotely located 
in a room overlooking the test sight. A description of the facility 
layout is given in Appendix B. 



4 
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V. CONCLUSIONS 



From the results of this investigation it is concluded that a 
time-shared digital-analog computer control system for aerodynamic 
heating is 

(1) capable of heat simulation control up to Mach 6. 

(2) silnple- 

(3) dependable. 

( 4 ) safe to operate. 

(5) relatively inexpensive. 
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Block Diagram of a Closed Loop Control System' for the 
Simulation of Kinetic Heating 

Fig. 1 
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Schematic Diagram of a Time Shaded' System for the 
Determination of "q^" 
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Reader* Read 




Example of multiplication of Vj_ by 11.7 



Reader switches- 8.0, 2.0, 1.0, .4, .2, and .1 are 
closed. Therefore; 

0 jR^ + R6 + R5 + ^3 R2 Rl) Vl 

or substituting values gives 



V 0 = -11.7 V, 



Diagram of Control System Multiplication Stage 

Fig. 3 



13 




O ^ £ 0 



Circuit Diagram of a Time Base Pulse- Generator 

Fig. 4 
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APPENDIX A 



Analog Summation 

The operational amplifier used, in the summation circuit has an input 
impedance which is assumed to he infinitely large. 

ft. 



V, 



R 1 

AAA- 



f 






-> V 



V, 



Since the input impedance is infinite, the current through terminal 
zero and currents tl 
equation may he written: 



A is zero and currents through R^ and R^ are equal. The following 



V - V V - V 
i a a o 



R, 



R, 



= 0 



( 1 ) 



also 



V = - M V 



( 2 ) 



Eliminating V from equations (l) and (2) gives; 



ft. 



( v i -Tf) ' - r i( v c + ^) 



(3) 



or rewriting 



V o E f M 

V. " R + R. ( 1 + MJ 
1 f l' ' 



(3a) 



If M is very large, the approximation is made: 



V - - 
o 



R f 

— v 

R 1 1 



(M 



For several inputs this can he extended. 
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-> V 




(5) 



Multiplication 

t is generated as a voltage signal by multiplying it by a reference 
voltage. 



^ — Output of -t (Step l) 



Ref • v 

-( 1 ) 

t is added to -t in a summer (See equation 5) to give (t - 1 ). 



aw 



Tape Reader 
Digital to 
Analog 



A 



,-t 



aw 



AW — i 



— VWVi 

-t> 



-> (t - t ) 
' aw s 



t >— / Vvv — 

s 



(Step 2) 



The analogous voltage, (t - t g ), is stored in a capacitor and 
appears as a continuous output of a peak follower. 



(t - t ) 
'aw s 





(Step 3) 



The output of the peak follower is fed back to the tape reader as a 
reference voltage where it multiplies h to give q . 



A-2 



(Step 4) 



J 



Tape Reader 
Digital to 
Analog 

A 



(t - t ) 
a v s 



— q = (t - t )h 
av s' 
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APPENDIX B: STANFORD UNIVERSITY RADIANT HEATING STRUCTURAL TEST FACILITY 



I. INTRODUCTION 

The need for aerodynamic heat simulation in structural testing of 
supersonic aircraft is veil recognized. Thus, the authors undertook 
the design and construction of a four channel radiant heating test facil- 
ity at Stanford University during the summer of 1962. The system vas 
designed to ieet the following requirements: a high intensity, rapid 

responso neat source; an adequate reflector system to direct the heat; 
an absorptive finish on the test specimen to permit maximum utilization 
of available heating power; a practical technique for surface temperature 
determination; a system of radiation measurement; and a computer control 
system to correctly vary the heat input to the test specimen. 

II . SYSTEM LOGIC 

Figure B-l is a block diagram of the control system for the radiant 
heating test facility. Voltages analogous to h and t are generated 
and then utilized by the computer along with a voltage analogous to t 
for solution of the Newtonian convective heat transfer equation, 
q = h(t Qw - tg) . Another voltage (q^) representing radiation from 
external sources (sun, atomic blast, etc.) is generated and then added 
to q in a summing circuit to give the desired heat flow signal, q^ . 
The« comparator output, c , is the difference between q^ and q^ , 
the measured heat flow. The control signal causes the power regulating 
device to vary the power supplied to the heaters. The heat flow measure- 
ment system closes the control loop, causing the correct heat input to 
the test structure to be maintained. The measurement of heat input to 
the test specimen is made by either of two systems. The operational 
details of these devices are explained briefly in Section III and in 
detail in References [l] and [ 2 ]. 

III. SYSTEM DETAILS 

Certain requirements are essential in a practical feedback control 
system for the simulation of aerodynamic heating. The requirements set 
for this facility are to : 
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(1) accurately compute desired heat flows for simulation of speeds 
up to Mach 6. 

(2) he capable of measuring the heat inputs to the test specimen. 

(3) be dependable. 

(4) be safe to operate. 

(5) be inexpensive as possible- 

(6) be simple to operate and program. 

The operations of the function generator and computer shown in 
Fig. B-l are combined in the digital to analogue computer control system- 
This unit consists of a punched tape reader, an operational amplifier 
and resistance bank, two peak followers and associated capacitors, a 
multi-bank stepping switch, a summer, and a chopper stabilized amplifier. 
(Fig. B-2) The details of the individual components are given 'in Ref- 
erence [l] and only the modus operandi is discussed here- 

Values of h and t are obtained from the Stefan-Boltzmann law. 

aw 

A binary coded program of these time dependent functions is made on a 
punched tape. This tape is placed on a reader and used in the time 
shared system as follows : an analogue voltage proportional to q^ is 

generated from a resistance network and a reference voltage in conjunction 
with an operational amplifier. This signal is temporarily stored. The 
next command is t which is treated in the same way. The output of 

the temperature transducers, t , is amplified by a chopper stabilized 

amplifier and then subtracted from the stored t analogue in a summer- 

The voltage difference (t - t g ) is used as the reference voltage for 

the next stage of the operation. This consists of giving the command, 
read h , and by virtue of the fact that the reference voltage is 
analogous to (t - t g ) the output of the system is analogous to the 
product h(t aw - tg) . Simultaneously, q and q_^ are added in a 

summer and this voltage, analogous to q, , is stored in a capacitor. 

Cl « 

These steps are repeated for channels two, three, and four. The digital 
to analog computer evaluates a new value of q^ for each channel every 
1.2 seconds. Thus, by time sharing, by digital to analog conversion, 
and by unique computer techniques, a simple and accurate means for 
obtaining a q^ history is realized. 
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There are two systems available for closing the control loop. Both 
systems measure the rate of heat flow to the test specimen and therefore 
may be used interchangeably. One system, the classic system, consists of 
special transducers. These essentially consist of thermocouples in con- 
junction with substantial heat sinks. In a radiant field the thermo- 
couples reach an equilibrium temperature proportional to the radiant 
intensity. Each transducer is two thermocouples connected so that their 
signals oppose. One is subjected to the radiation from the heaters and 
the other to the reflected radiation. The output is therefore proportional 
to the difference between the incident and reflected radiation. This is 
exclusive of convective losses. The heat sinks for the transducers 
are cooled metallic blocks. The instruments are attached directly to the 
sides of the heater arrays. They are calibrated in a calorimetric manner. 

The other system, one especially devised for this project, consists 
of eighteen thermocouples arranged to obtain temperature profiles in the 
test specimen. A high speed stepping switch sequentially connects the 
thermocouple signals to an analogue computer which solves the appropriate 
heat flow equation. The computer is time shared by the 4 stations. The 
latter system is more versatile than the former -since it measures net 
heat flow irrespective of the nature of the source. Thus convective losses 
in simulation are accounted for. An additional refinement to the computer 
permits compensation for radiant losses which the structure would normally 
expedience in flight. Figure B-3 is a photograph of the feedback com- 
puter system. 

Voltages in the control system are scaled such that at any time the 
prescribed q^ level is twice the anticipated q^ level. A summing 
circuit subtracts q^ from q^ and the output (c) is used to control 
the power regulating devices. Due to the two-to-one scaling the control 
signal consists of a set-point voltage plus an error voltage. Thus, the 
heaters are controlled in such a manner that the desired and actual heat 
inputs are equalized. This occurs when the error voltage level converges 
to zero and the control signal converges to the set-point voltage. 

Main power is supplied through two 300 KVA, 480 volt, 3 phase, 

General Electric transformers. The power regulating devices consist of 
75 KVA, 480 volt, 3 phase General Electric saturable core reactors and 
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their associated control units. Transient operation of the power equip- 
ment on an overload duty cycle makes possible peak outputs of two to 
three times rated value- Heating capacity of more than one megawatt is 
therefore feasible. Figure B-4 is a photo of the saturable core reactor 
units . 

The basic heating elements are General Electric T-3 infra-red lamps. 
(See Fig. B-5 ) They consist of a tungsten filament supported by tanta- 
lum disc spacers in a sealed quartz tube 3/8 of an inch in diameter and 
12 inches in length. Each unit is rated at 1000 watts at 220 volts. 

They are capable of withstanding the k80 volt peak output of the power 
supply and under these conditions develop approximately 3 kilowatts. 

The heaters are arranged in groups of six. Figure B-6 is a photo- 
graph of one such unit. The heater arrays are connected in delta to the 
three phase power from the saturable core reactors. Each unit covers 
about one third of a square foot. Thus, the maximum radiant intensity 
possible with the present arrangement is about 60 kilowatts per square 
foot- Two thousand watt heaters of the same external size are avail-' 
able when higher intensities are required. Peak output of these is 6000 
watts per heater. The heater arrays are quite adaptable and may be 
bolted to any suitable frame to surround a test structure. Special 
arrangements may easily be devised for special purposes, such as the 
heating of curved surfaces. 

F-igure B-7 shows an overall view of the test facility. The main 
items of equipment and the channels of information flow are outlined on 
the photograph. All control, computing and recording equipment is located 
in the control room. This room is on the second floor of the laboratory 
overlooking the test site. Large windows provide for visual observation 
of the structure under test. The remote location of the control room 
was dictated by safety considerations and provides a convenient, compact 
center for operation of the entire facility. 

IV. CONCLUDING REMARKS 

The design, construction, and operation of this test facility has 
proved that the concept of component time-sharing is practical in certain 
types of control systems. It is concluded that such time-sharing is an 
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economical proposition in which savings in equipment costs are balanced 
against slightly reduced system performance. The facility is such that 
it can be readily installed in any small research laboratory where cost 
is of prime importance. 
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Bloclr Diagraxa of a Closed Loop Control System- for tho 
Simulation of Aerodynamic Heating 



Fig. B-l 
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Digital to Analog Computer Control Feedback Computer 

Fig. B-2 Fig. B-3 




Saturable Core Reactors 
Fig. B-4 




Infra-Red Heater 
Fig. B-5 



Heater Array 
Fig. B-6 
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Overall View of Test Facility 
Fig. B-7 




Control Room 
Fig. B-8 
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APPENDIX C 



TEST DATA 
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TABLE I 

DUPLICATION AND ACCURACY OF MULTIPLICATION 



Reference Voltage 


Multiplier 


Actual Prod. 


°jo Error 


1.50 


volts 


5.60 


5.35 


volts 


0.9 


1.88 


volts 


3.60 


6.70 volts 


1.0 


2.02 


volts 


3.10 


6.26 volts 


0.0 


3.10 


volts 


3.10 


9.60 volts 


0.1 


O 

to 


volts 


3.00 


12.7 


volts 


0.9 


6. 20 


volts 


3.00 


18.3 


volts 


1.6 


8.42 


volts 


2-95 


24.4 


volts 


1.6 


10.4 


volts 


2.95 


30.3 


volts 


1.3 


12 ’.4 


volts 


2.60 


32.0 


volts 


0.6 


13.8 


volts 


2.60 


35.3 


volts 


1.4 


15.2 


volts 


2.20 


33.8 


volts 


1.2 


19.0 


volts 


2.20 


42.1 


volts 


0.5 




Duplication of Above 


Test after 15 Minutes : 




1.50 


volts 


3.60 


5.35 volts 


0.9 


1.88 


volts 


3.60 


6.70 volts 


1.0 


2.02 


volts 


3.10 


6.28 volts 


0.3 


3.10 


volts 


3.10 


9.61 volts 


0.0 


' 4.30 


volts 


3.00 


12.8 


volts 


0.8 


, 6.20 


volts 


3.00 


18.4 


volts 


1.0 


Oo 

V- 

ro 


volts 


2.95 


24.5 


volts 


1.0 


10.4 


volts 


2.95 


30.5 


volts 


0.6 


12.4 


volts 


2.60 


31.9 


volts 


1.0 


13.8 


volts 


2.60 


35.5 


volts 


0.8 


15.2 


volts 


2.20 


33.5 


volts 


0.3 


19.0 


volts 


2.20 


42.5 


volts 


1.4 
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0 2 A 6 8 10 12 14 16 18 20 22 24 26 28 30 32 



•rims (t 30 Minutes) seconds > 



A s trip- chart recording of an arbitrary function 
analogous to q d starting at time t and at time t -1-30 minutes 



Fig. C-l 
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Tine t (.5 era /see) 




Tirr.o t (.3 cm/sec) 



3tr' p- chart recording of constant o , generated on each 

of four channels 

Fig. C-2 
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